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ABSTRACT 

A major recent development in observational cosmology has been an accurate mea- 
surement of the luminosity distance-redshift relation out to redshifts z = 0.8 from 
Type la supernova standard candles. The results have been argued as evidence for 
cosmic acceleration. It is well known that this assertion depends on the assumption 
that we know the equation of state for all mass-energy other than normal pressureless 
matter; popular models are based on either the cosmological constant or on the more 
general quintessence formulation. But this assertion also depends on a number of other 
assumptions, implicit in the derivation of the standard cosmological field equations: 
large-scale isotropy and homogeneity, the flatness of the Universe, and the validity of 
general relativity on cosmological scales (where it has not been tested). A detailed ex- 
amination of the effects of these assumptions on the interplay between the luminosity 
distance-redshift relation and the acceleration of the Universe is not possible unless 
one can define the precise nature of the failure of any particular assumption. However 
a simple quantitative investigation is possible and reveals a number of considerations 
about the relative importance of the different assumptions. In this paper we present 
such an investigation. We find that the relationship between the distant-redshift re- 
lation and the sign of the deceleration parameter is fairly robust and is unaffected if 
only one of the assumptions that we investigate is invalid so long as the deceleration 
parameter is not close to zero (it would not be close to zero in the currently-favored 
r^A = 1 — r^mattcr = 0.7 Or 0.8 Uuivcrse, for example). Failures of two or more assump- 
tions in concordance may have stronger effects. 
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1 INTRODUCTION AND BACKGROUND 



One of the most significant cosmological results of recent 
times has been an accurate determination of the luminosity 
distance - redshift relation (Perlmutter et al. 1999, Schmidt 
et al. 1998, Riess 1999, Zehavi & Dekel 1999, Riess et 
al. 2000, Burrows 2000), from measurements of Type la su- 
pernova standard candles out to redshifts z ~ 0.8. The mea- 
surements are highly inconsistent with a flat purely matter- 
dominated Universe and have consequently been regarded 
as evidence for cosmic acceleration. 

But this assertion depends intricately on other assump- 
tions made, and does not necessarily follow from the mea- 
surement of the luminosity distance-redshift relation. Sup- 
pose that we can describe the expansion of the Universe 
by a scale-factor R{t) which is a function of cosmic time t. 
The redshift z — R{to) / R{t) — 1 is then also a function of t 
(here to is the current age of the Universe). The luminosity 
distance-redshift relation is (e.g. Weinberg 1972) 



dL{z) = R{to){l + z)Sk 



cdR 

dR n 



(1) 



where c is the speed of light, k is an additional parameter 
describing the curvature of the Universe (fc = -1-1 for a closed 
Universe, k = for a flat Universe, and k — —1 for an 
open Universe), and S is defined so that S'-|-i(a;) = sinx, 
So{x) — X, and S-i{x) — sinha;. The deceleration parameter 
is (e.g. Weinberg 1972) 



go = 



(2) 



If the universe is accelerating qo < 0. It foUows from these 
equations that a function R{t) constrained from a measure- 
ment of dL{z) does not uniquely define a value of qo. In fact, 
it does not even uniquely define the sign of qo ■ This happens 
because go depends on the second time derivative of R{t) but 
dhiz) only depends on lower order derivatives and their in- 
tegral. As a simple (albeit unphysical) example, consider a 
function R(t) that can be expressed as a polynomial in t 



© 0000 RAS 



2 Neil Trentham 



with coefficients a„ such that R{t) = Sa„t" where the nth 
term contributes significantly more than the (n + l)th one. 
If one then defines another function R!{t) that is identical 
to R{t) except that the sign of the 02 coefficient is opposite, 
then both R{t) and R'(t) generate almost identical dL(z) 
relations but deceleration parameters with opposite signs. 

It is conventional to draw a connection between the 
dhiz) relation and the sign of the deceleration parameter as 
follows. One assumes the Friedmann-Robertson- Walker met- 
ric (which follows from the requirement of isotropy and ho- 
mogeneity on the largest scales), Einstein's theory of general 
relativity, zero spatial curvature (fc = 0) and an equation 
of state so that the only contributions to the stress-energy 
tensor in general relativity comes from normal matter and 
the vacuum energy (the cosmological constant). This latter 
assumption is frequently examined (Caldwell, Dave & Stein- 
hardt 1998; Garnavich et al. 1998; Efstathiou 1999; Wang et 
al. 2000; Maor, Brustein & Stoinhardt 2000) since it has no 
theoretical basis. If those assumptions are made, the model 
forces a particular form of R{t) and consequently a par- 
ticulax form of cLl^z) and go in terms of Q,m and Q,a, the 
cosmological densities of matter and the vacuum. The de- 
celeration parameter here has the particularly simple form 
go = f2M/2 — f^A. Given the constraints imposed by the 
supernova measurements, it is found that the only pairs of 
Q,M and Q,k that are allowed force go < 0. Therefore the 
Universe is argued to be accelerating. 

But any of the assumptions may turn out to be invalid. 
Aside from the much-examined assumptions regarding the 
equation of state turning out to be untrue, it is also possible 
that fc 0, that very-large-scale inhomogeneitios may exist 
(Celerier 2000), or even that general relativity may turn out 
to be invalid on the very largest scales (where is has not been 
directly tested). If any or all of these turn out to be true, 
we might no longer require an accelerating Universe given 
the current (or even a far more precise) measurement of the 
dhiz) relation; this possibility is investigated in a simple 
quantitative way in this paper. 



2 FRIEDMANN MODELS 

The formulation that is normally used to define the cosmo- 
logical parameters can be conceptually (albeit arbitrarily) 
divided into the following four steps. 

(i) The cosmic geometry is described by the four- 
dimensional Riemann curvature tensor. If isotropy and ho- 
mogeneity are assumed, the resulting constraints on the Rie- 
mann tensor force the space-time metric to have a specific 
form, given by the Priedmann-Roberston- Walker metric. Im- 
plicit in this form is the parameter k, which takes the value 
or ±1. The existence of this parameter k follows from 
the assumption that the Ricci scalar (the one-dimensional 
fully-contracted Riemann tensor) is not a function of either 
time or of position, which is required by homogeneity. If this 
assumption about homogeneity and isotropy is made, then 
observations of the position /peak of the first Doppler peak 
in the microwave background spectrum (e.g. de Bernardis 
et al. 2000) that suggest a value of the total cosmological 
density Otot close to 1 in units of the critical density then 
suggest that A; = 0. 

(ii) Prom the metric, one can then compute explicitly the 



curvature tensor and its contractions. From the equations 
of general relativity, one can then compute the cosmologi- 
cal field equations. Note that general relativity has not been 
verified observationally to high precision on these cosmolog- 
ical scales (although it is clearly very successful on smaller 
scales - see Section 8 of Weinberg 1972) so its application 
here is very much an assumption. 

(iii) if we then adopt an equation of state that relates den- 
sity to pressure for every contribution of mass density in 
the field equations, we can then derive R(t) as a function of 
these contributions. Some common strategies are to assume 
that (a) the only contribution comes from normal pressure- 
less matter; (b) part of the contribution comes from normal 
matter and all the remaining part comes from a cosmologi- 
cal constant; (c) part of the contribution comes from normal 
matter and all the remaining part comes from some mate- 
rial whose equation of state is defined by the quintessence 
formulation of Caldwell et al. (1998). 

(iv) from measurements of luminosity distance-redshift re- 
lation we constrain combinations of parameters that arise 
from the considerations in (iii) above. For example, if one has 
specified exactly two forms of mass density, both of which 
have known equations of state, then one can constrain the 
relative mass densities of the two components (this would be 
true in case (b) above) . If one makes a further assumption, 
like requiring a flat Universe (k = Q), then one has another 
constraint and in conjunction with the luminosity distance- 
redshift relation can provide a precise determination of the 
model parameters. Alternatively any related measurement 
like the value of Zpeak will have a similar effect (recall from 
(i) that in the context of these models, the value of /peak that 
was measured is operationally similar to requiring A; = 0). 
In practice, the constraints on the model parameters in case 
(b) are tight since the constraints on the mass densities in 
normal matter and in the cosmological constant are almost 
orthogonal (see Efstathiou et al. 1999). From the permitted 
values of the contributions of these mass densities, we can 
then compute R{t) and hence go, within the errors provided 
by the measurements. 



3 MODIFICATIONS TO FRIEDMANN 
MODELS 

Clearly there are lots of assumptions involved. It is there- 
fore useful to examine each in detail and to see how re- 
laxing each affects the interplay between c?l(z) and go. As- 
sessing the general case is mathematically complicated and 
unconstrained. However, it is possible to compute the con- 
sequences of relaxing the assumptions in simple well-defined 
ways, and this is done in Table 1. The approach is as follows: 

(1) One of the assumptions described in the previous sec- 
tion is relaxed (given in the first column of Table 1) and a 
modification of the relevant equation is assumed. Each ad- 
justment is described in terms of a parameter e, which is a 
measure of how big the modification is. 

(2) The relevant modified Friedmann equation is derived 
(the second column of Table 1). 

(3) From this equation, we then compute the value of rfi (2: = 
1) for various fiducial models, listed in Table 2 (note that 
Model C is consistent with the observations described in 
Section 1, but Models A and B are not). These represent 
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Table 1. Simple modifications to the Priedmann models 



Modification and form 



New Priedmann equation" 



Isotropy and homogeneity 

TZ ~ sin^ e 
TZ ~ sin^ d> 



as standard 

as standard 

(-e - e cos 26» + 2 sin^ 6») e sin-2+^ e + 
4r2(l + 3sin^ 9 - sin^^ 6>)i?2 = 32TrGpr^R'^ 

(2 — e - ecos 24>) e coscc^e sin^^+^ (p + 
4r2(l + Ssin-^ (j> - sin^^ </>)i?2 = 327rGpr^R^ 



Density (curvature) 
+ Ha = 1 + e 



as standard 



General relativity 



R \-3(l+e) 



Equation of state 

P = tp (quintessence'^) H'^ = ^mH^ {^) '^ + (^- ^m)HI [^) 



The symbols have the following meanings: e is a parameter, defined individually in each case, that describes how wc arc modifying 
the model; {t,r,6,(j>) arc the time, radial, and two angular coordinates used to define the space-time metric g^,^; the two-dimensional 
contraction of the full four-dimensional Ricmann curvature tensor is the Ricci tensor R/j^v and its contraction to one dimension is the 
Ricci scalar TZ; T^i^ is the Einstein stress-energy tensor, implicit in which are the density p and pressure P; R is the cosmic scale factor, 
which depends on t, R is its time derivative and = ilo and Hq (the Hubble constant) are the present-day values of R and H; Qm 
is the cosmological density of normal matter obeying the equation of state P = in units of the critical density; Qa = ^/■^Hq is the 
cosmological density in the cosmological constant; G is the gravitational constant. Other symbols are defined elsewhere in the text; the 
speed of light c is set to unity in all the equations in the table. 

Assuming fe = and A = 0. 

^ In order to fully define a quintessence model we need to state both the equation of state of the quintcsscnt material (parameterized by 
e here) and state the density of this material. For the models wc consider here, wc fix the density of normal matter to be Qm = 0.2 in 
units of the critical density and the density of the quintcsscnt material to be Sle = 0.8 so that the total density is f^TOT = 1. The value 
of Qm = 0.2 is suggested by dividing the value of Hbaryon derived from big bang nucleosynthesis constraints (e.g. Smith, Kawano & 
Malaney 1993) by the baryon fraction in rich galaxy clusters like Coma (White et al. 1993). It is also the value suggested for ildark matter 
by (see Trimble 1987) assigning an amount of dark matter to each luminous galaxy (e.g. using the correlations of Kormendy 1990) and 
integrating over the galaxy luminosity function (e.g. Ellis et al. 1996). Neither of these methods depend on the presence or absence of a 
form of mass density that does not behave according to the equation of state P = 0, like a cosmological constant. The value of S^xoT = 1 
is suggested by the measurements of the cosmic microwave ba<;kground by de Bernardis et al. (2000). 



Table 2. Fiducial models 



Model 


Parameters 


k 


dL{z = l)/cH-'' 


90 


A 


CIm = 1 





1.17 


0.50 


B 


Om = 0.2 


-1 


1.41 


0.10 


C 


CIm = 0.2, Ha = 0.8 





1.65 


-0.70 
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Table 3. Eflfects of modifications on observables 



Modification and form 



Values of e that 
are consistent with a 
measurement of di^{z = 1) 
that is accurate to 10% 



Values of e that 
are consistent with the 
fiducial value of qo 
to within 10% 



Values of e that 
keep unchanged 
the sign of qo 



Isotropy and homogeneity 

-R. ~ sin« 9 
TZ ~ sin^ d> 



any e 
any e 
depends on (r, 0) 
depends on (r, 0, cp) 



any e 
any e 
depends on (r, 6) 
depends on (r, 0, </>) 



any e 
any e 
depends on (r, 6) 
depends on (r, 0, ^) 



Density (curvature) 



1 + e 


A: 


-0.45 < e < 0.63 


-0.10 


< 


e 


< 0.10 


e > -1.00 




B: 


-0.30 < e < 0.41 


-0.02 


< 


e 


< 0.02 


e > -0.20 




Cl^ 


-0.18 < e < 0.54 


-0.14 


< 


e 


< 0.14 


e < 1.40 






-0.50 < e < 0.18 


-0.07 


< 


e 


< 0.07 


e > -0.70 



General relativity 



87rGr^.(r«)^ 


A: 


-0.19 < e < 0.22 


-0.03 


< 


€ < 0.03 


e > -0.33 




B: 


-0.89 < e < 0.52 


-0.03 


< 


e < 0.03 


e > -0.33 




C: 


-0.39 < e < 0.36 


-0.23 


< 


e < 0.23 


e < 2.33 



Equation of state 

P = tp^ A: -0.295 < e < 0.27 -0.04 < e < 0.04 e > -0.42 

B: -0.79 < e < -0.21 -0.34 < e < -0.33 e > -0.42 

C; -1.41 < e < -0.65 -1.06 < e < -0.94 e < -0.42 



* For the fiducial model C, modifications that are introduced can be considered as modifications to the behaviour of the matter component 

(line CI) or the cosmological constant component (line C2). Results are presented for both. 

^ Note that the value of e = is not "permitted" for models B and C. This is because in neither case can the fiducial model be described 
as an unmodified version of the model under consideration here. In fact, fiducial model C is exactly equivalent to a modified model with 
e = —1; this value of e is of course at the centre of the permitted range. The reason that models A, B, and C give different constraints 
on e in the second and third columns is because differrent fiducial models have been adopted. In the fourth column, where the limiting 
e does not depend on the fiducial model, this limiting value is the same in all three cases. 



the case e = 0. We then ask: what values of e generate value 
of di(« = 1) that differs by less than ten per cent from 
this fiducial value? These values arc presented in Table 3. 
A measurement of the luminosity distance at z = 1 that 
is precise to 10% requires observations slightly deeper and 
more precise than those of Perlmutter ot al. 1999, but should 
be attainable with larger samples of supernovae observed 
with large telescopes in the near future. 

(4) We also compute the values of e that leave the sign of the 
deceleration parameter within 10 % of its fiducial value and 
that leave the sign of the deceleration parameter unchanged 
(the third and fourth columns of Table 3). 

If the values of e allowed in the second column of Table 3 
all satisfy the condition given in the fourth column of Table 
1, then we can conclude that a measurement of dL{z = 1) 
precise to 10% does indeed constrain the sign of go (for ex- 
ample, this is true of the case labeled "General Relativity" ) 
for the relevant fiducial model. 



An example of this calculation follows. Consider the last 
row in the "density (curvature)" section of the table i.e. the 
line labelled "C2" . Given the Priedmann equation, equation 

(1) now becomes 



where k ■ 
and 



e\Ho 

+1 if e > 0, A; 



-f 

Jo 



■1 if € < and k 

dz' 



(3) 

if e = 0, 
(4) 



^(1 + z')^{l + 0.2z') + (0.8 + e)z'{2 + z') ' 

If e = (fiducial model C), then dz,(« = 1) = 1.65 c/Hq. 
For a value within 10 per cent of this (IASc/Hq < dL{z = 
1) < 1.81c/i^o), Equation (3) requires -0.50 < e < 0.18. 
Equation (2) now becomes 

qo = -0.7 - €. (5) 
For the fiducial model qo = —0.7 and is negative. Values of 



© 0000 RAS, MNRAS 000, 000-000 



Distance measurements as a probe of cosmic acceleration 5 



e that result in a value of go within 10% of this value are 
—0.07 < e < 0.07, from Equation (5). To keep qo negative, 
Equation (5) requires e > —0.70. Therefore any changes in e 
that are consistent with a measurement of dL{z = 1) precise 
to 10 per cent do not change the sign of qo, at least for this 
fiducial model and form of perturbation. 

Such numbers, although interesting, are of limited use 
since the form of the modification that we use is so specific. 
We would, however, make the following two comments: 

(i) the choice of the fiducial model is important. For fiducial 
models like Model B which have qo close to zero, a small 
change in dL{z = 1) could be consistent with a model in 
which the sign of qo changes. However, for model C, which 
is consistent with current observation, the fiducial value of 
go is far enough away from zero that all of the modifica- 
tions considered here do not change the sign of go without 
changing dL{z = 1) at the 10% level. 

(ii) small dependences of the Ricci scalar on any metric co- 
ordinates introduce small changes in the Friedmann equa- 
tion. Larger dependances that are power laws in either time 
or in the radial position coordinate also have negligible ef- 
fect. Larger dependencies having other functional forms will 
have a bigger effect but these will probably need to be fairly 
finely tuned if they are going to affect the relationship be- 
tween dL{z) and go and not have any appreciable signature 
in either the galaxy distribution on large scales or the cosmic 
microwave background; 

(iii) invalidating general relativity in the way defined in Ta- 
ble 1 does not seem to have a strong effect on the inter- 
play between dL{z — 1) and go. Values of e large enough to 
change the sign of go would cause a very substantial change 
in dL{z = 1) (except for Model B, but see point (i)); 

(iv) Modifying the equation of state according to the 
quintessence formulation has some effect on the interplay 
between dL{z — 1) and go. Maor et al. (2000) consider the 
case of a time-varying equation of state, which is outside 
the bounds of the simple mathematical treatment presented 
here, and find that relaxing the assumption of a constant 
equatio n of state can have very profound implications for 



using dc(z) to determine go and more generally the fate of 
the Universe. 

In these calculations, we have concentrated on measure- 
ments of dL{z) at z = 1. In principle, measurements at 
higher redshifts would help since ^ g^^"^ increases as z in- 
creases (formally, this is true only if the underlying cosmol- 
ogy is known a priori; see also the discussion on this point 
in Maor et al. 2000). But one would then have additional 
system atic effects that would need to be understood, like 



K-corre ctions of Type la supernovae, complications due to 



dust extinction (Riess et al. 2000) and possibly even gravi- 
tational lensing by matter along our line of sight (Metcalf & 
Silk 1999). More precise measurements at z = 1 and lower 
redshifts would also help, although these may be fundamen- 
tally limited by intrinsic scatter in the peak luminosity of 
Type la supernovae (see the references in the first paragraph 
of Section 1). 



the value of go is close to zero. For the SIa = 1 — JIm = 0.8 
model currently favoured by observation, the value of go is 
far enough away from zero that this concern is not valid. 
This analysis is, however, limited, in that it only consid- 
ers modifications to the Friedmann models of a very spe- 
cific type, and those individually. The joint effects of two 
or more modifications might be more dramatic, but at this 
stage there are too many unconstrained ways to formulate 
these effects for a detailed analysis to be productive. In 
the long term, measurements like those of the cosmic mi- 
crowave background (Efstathiou et al. 1999), gravitational 
lensing statistics (e.g. Waga & Frieman 2000) and galaxy 
number counts (Newman & Davis 2000) will provide other 
constraints. 

In order to address formally the question as to whether 
or not the Universe is accelerating, it will be necessary to 
have measurements which probe the second time derivative 
of R{t). However other, less direct, measurements like the 
ones listed in the previous paragraph and the supernova 
measurements will provide powerful hints if the association 
between the luminosity distance-redshift relation and the 
sign of go really does turn out to be only weakly model- 
dependant. 
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4 SUMMARY 

Therefore the relationship between dL{z = 1) and the sign 
of go is fairly robust given the analysis presented here unless 
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